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Quantum Finance represents the synthesis of the techniques of quantum theory 
(quantum mechanics and quantum field theory) to theoretical and applied finance. 
After a brief overview of the connection between these fields, we illustrate some 
of the methods of lattice simulations of path integrals for the pricing of options. 
The ideas are sketched out for simple models, such as the Black-Scholes model, 
where analytical and numerical results are compared. Application of the method 
to nonlinear systems is also briefly overviewed. More general models, for exotic or 
path-dependent options are discussed. 



1. Introduction 

^ Financial markets have undergone a tremendous growtli in the last 
decades and in order to meet the need of customers new and complex 
financial instruments have been developed. Risk assessment models and 
the quantification of returns, given the huge amount of trading involved 
worldwide, requires more sophisticated approaches than in the past. Non- 
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linearities play a key role in all of this, and, from this side, the field is 
probably largely unexplored. 

In general, intensive numerical simulations and fast algorithms are 
needed to obtain useful results. Analytical results are limited, except for 
simple models such as the Black-Sholes model and other similar models. 

Use of a path integral formulation has some advantages. First, it is in 
close relation to the lagrangean description of diffusion processes, second, 
it opens the way to the use of quantum mechanical methods. 

We will introduce the field at a non-expert level through a crash course 
(next few sections). Then we will come to brefly illustrate where nonlinear- 
ities appear and review some of the simplest equations one can write down, 
generalizing the Black-Scholes model. We will then proceed to discuss the 
path integral formulation of the Black-Scholes and the model of barrier 
options, resorting to a lagrangian path integral formulation. Strategies to 
solve the path integral are then briefly presented, together with some re- 
sults. Detailed simulations, algorithms of analysis and related applications 
will be presented elsewhere ^ . The review sections are standard knowledge 
in the field, and are based on ^. 

2. A view on Theoretical Finance 

2.1. Simulating the Complexity of Finance 

The simulation of financial markets can be modeled, from a theoretical 
viewpoint, according to two separate approaches: a bottom up approach 
and (or) a top down approach. 

For instance, the modeling of financial markets starting from diffusion 
equations and adding a noise term to the evolution of a function of a 
stochastic variable is a top down approach. This type of description is, 
effectively, a statistical one. 

A bottom up approach, instead, is the modeling of artificial markets 
using complex data structures (agent based simulations) using general up- 
dating rules to describe the collective state of the market. The number 
of procedures implemented in the simulations can be quite large, although 
the computational cost of the simulation becomes forbidding as the size 
of each agent increases. Readers familiar with Sugarscape Models and the 
computational strategies based on Growing of Artificial Societies ^ have 
probably an idea of the enormous potentialities of the field. However, one 
would expect that the bottom up description should become comparable to 
the top down description for a very large number of simulated agents. 

The bottom up approach should also provide a better description of ex- 
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treme events, such as crashes, collectively conditioned behaviour and mar- 
ket incompleteness, this approach being of purely algorithmic nature. A 
top down approach is, therefore, a model of reduced complexity and fol- 
lows a statistical description of the dynamics of complex systems (for an 
introduction see 

2.2. Forward, Futures Contracts and Options 

Let the price at time t of a security be S{t). A specific good can be traded 
at time t at the price S(t) between a buyer and a seller. The seller (short 
position) agrees to sell the goods to the buyer (long position) at some time 
T in the future at a price F(t,T) (the contract price). Notice that contract 
prices have a 2-time dependence (actual time t and maturity time T). Their 
difference t = T — t is usually called time to m,aturity. Equivalently, the 
actual price of the contract is determined by the prevailing actual prices 
and interest rates and by the time to maturity. 

Entering into a forward contract requires no money, and the value of the 
contract for long position holders and strong position holders at maturity 
T will be 

{-ir{S{T)-F{t,T)) (1) 

where p = for long positions and p = 1 for short positions. Futures 

Contracts are similar, except that the after the contract is entered, any 
changes in the market value of the contract are settled by the parties. 
Hence, the cashflows occur all the way to expiry unlike in the case of the 
forward where only one cashflow occurs. They are also highly regulated 
and involve a third party (a clearing house). Forward, futures contracts 
and, as we will see, options go under the name of derivative products, since 
their contract price F(t, T) depend on the value of the underlying security 
S{T). 

To complete this crash course on financial instruments we need to define 
options. Options are derivatives that can be written on any security and 
have a more complicated payoff function than the futures or forwards. For 

example, a call option gives the biiyer (long position) the right (hut not the 
obligation) to buy or sell the security at some predetermined strike-price at 
maturity. A payoff function is the precise form of the price. Path dependent 
options are derivative products whose value depends on the actual path 
followed by the underlying security up to maturity. In the case of path- 
dependent options, since the payoff may not be directly linked to an explicit 
right, they must be settled by cash. This is sometimes true for futures and 
plain options as well as this is more eflacient. 
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3. Langevin Evolution 

In the top down description of theoretical finance, a security S{t) follows a 
random walk described by a Ito-Weiner process (or Langevin equation) as 

= (l,dt + aR{t)dt, (2) 

where R(t) is a Gaussian white noise with zero mean and uncorrelated 

values at time t and t' {R{t)R{t')) = S{t — f/). </) is the drift term or expected 
return, while a is a constant factor multiplying the random source R{t), 
termed volatility. 

As a consequence of Ito calculus, differentials of functions of random 

variables, say f{S,t), do not satisfy Leibnitz's rule, and for a Ito-Weiner 
process with drift (2) one easily obtains for the time derivative of f{S,t) 

df df 1 2n2d^f ,:odf , adf 

di=dt+r^dS-^+'^^dS^''^ds'^- 

The Black-Scholes model is obtained by removing the randomness of the 
stochastic process shown above by introducing a random process correlated 
to (3). This operation, termed hedging, allows to remove the dependence on 
the white noise function R{t), by constructing a portfolio 11, whose evolution 
is given by the short-term risk free interest rate r 

dU , , 

-dt='''- 

A possibility is to choose 11 = / — ^S". This is a portfolio in which the 
investor holds an option / and short sells ^ an amount of the underlying 
security S proportional to A combination of (3) and (4) yields the 
Black-Scholes equation 

There are some assumptions underlying this result. We have assumed ab- 
sence of arbitrage, constant spot rate r, continuous balance of the portfolio, 
no transaction costs and infinite divisibility of the stock. 

The quantum mechanical version of this equation is obtained by a 
change of variable S = e^, with x a real variable. This yields 

%=H.sf (6) 



''short selling of the stock should be possible 
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with an Hamiltonian Hbs given by 

Notice that one can introduce a quantum mechanical formahsm and inter- 
pret the option price as a ket |/) in the basis of |a;), the underlying security 
price. Using Dirac notation, wc can formally reinterpret f{x,t) = {x\f{t)), 
as a projection of an abstract quantum state \ f{t)) on the chosen basis. 

In this notation, the evolution of the option price can be formally written 
as I/, t) = e*^|/, 0), for an appropriate Hamiltonian H. 

In the presence of a stochastic volatility, the description is more involved, 
but also more interesting. 

In general, the description of these processes is driven by two correlated 
white noise functions R\ and R2 

^ = A + Ml^ + CV^"-Ri 
at 

^=(i)S + a^ + liV + CV''R2 (8) 

with V = ^fa and {Ri{t)R2{t')) = 1/p S{t — f/), p being the correlation pa- 
rameter. However, since volatility is not traded in the market (the market 
is said to be incomplete), perfect hedging is not possible, and an additional 
term, the market price of volatility risk /3{S,V,t,r), is in this case intro- 
duced. f3 can be modeled appropriately. In some models ^, a redefinition 
of the drift term /u in (8) in the evolution of the volatility is sufficient to 
hedge such more complex portfolios, which amounts to an implicit choice of 
f3{S, V, t, r). We just quote the result for the evolution of an option price in 
the presence of stochastic volatility, which, in the Hamiltonian formulation 
are given by 

%=HMGf (9) 



dt 



with 



Hmg- \t ^Jg^ ^Xe +^i j ^ 

which is nonlinear in the variables x = log(5) and y = log(y). For general 
values of the parameters, the best way to obtain the pricing of the options 
in this model is by a simulation of the path integral. 
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4. Monte Carlo Simulations of option pricing 

Simulations of the functional integral are rather straightforward and we 
should omit any detail about them since they have been known ever since 
by the high energy physics community. However, to reach out to a less 
specialized audience, we will provide a simple illustration of the method. 

Once the model is given, one determines the underlying lagrangean. We 
assume a discretization of the time to maturity r in intervals e = t/N, with 
N an arbitrary (large) integer. 

For instance, for the Black Scholes model one gets the action 

JV 

SBS = eJ2LBs{i) (11) 



with 
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where we have introduced discretized positions (a^j) for the variable x = 
log(S') which identifies the quantum mechanical state of the system. We 
will refer to it as to the stock price. The propagator for the stock price is 
given by the pricing kernel 

Pbs{x,x',t) = J DXBse'^^' 
= (a;|e-^^«^|a:') 

(13) 

with 

/ DXbs = n[=o r dx{t). (14) 

J J — CO 

For barrier options there is an analogous procedure, except that now we 
need to introduce a generic potential V{x) in the corresponding Hamilto- 
nian 

The pricing kernel is the fundamental quantity to compute using the func- 
tional integral. Related attempts can be found in the literature 

For this purposes, we have used a standard Metropolis algorithm. If 
thermalization is slow, it is possible to resort to use sequentially Metropolis 
updates and cluster updates. The latter is an update for the embedded Ising 
dynamics in the lattice variables Xi/\xi\ (Swendsen-Wang, Wolff), and is 
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included in for a faster generation of the thermalized paths of the stock 
price x{t). 

For processes involving a stochastic volatility {y = log(V)) the expres- 
sion of the path integral is more complicated and can be found in ^. From 
now on we will just consider the case of a constant volatility. 

If we denote by g{x, K) the payoff function, with a strike price K, in 
this case the value of the option (its price) is given by the Feynman-Kac 
formula 

/oo 
dx' {x\e-^^-*^"''^\x')g{x' ,K). (16) 
-oo 



In actual simulations, it is convenient to compute directly the option price 
rather than the propagator itself. The simulation is done by taking the 
initial point x fixed, and letting the final point evolve according to its 
quantum dynamics. In this way a path {x,x') is generated. After the first 
thcrmalization, x' is allowed to undergo quantum fluctuations, at fixed 
X. Each x' is then convoluted with the payoff function and an average is 
performed. Finally, this procedure is repeated for several x values, so to 
obtaint the option price at time to maturity r. 

Figs. 1, 2 and 3, illustrate some simple results obtained by the monte 
carlo mehtod. For illustrative purposes, we show the behaviour of the 
Black-Scholes model. Fig. 1 shows a typical thermalized path, generated 
from a given initial value x (at current time t = t) assuming a maturity of 
300 days, while in Fig. 2 we have plotted several path for different starting 
values X of the stock at current time r. We have chosen an interest rate 
r = 0.05 and a 12 percent volatility a. Finally, in Fig. 3 we compare the 
analytical and the numerical evaluation of the Black-Scholes option price 
with a low resolution for (16), in order to separate the two curves, which 
otherwise would overlap completely, in order to illustrate the convergence 
of the Metropolis algorithm. 

Barrier options can be analized similarly, equivalently, by this method 
or by the Langevin method. We show in Fig. 4 the evaluation of the price 
of the option using the Langevin method in the presence of a step potential 
sitting at a value of the stock price given by 2:0 = log(S'o), with 5*0 = 100. 
Compared to the Black-Scholes now the price has been discounted. 

Applications of the method to the determination of various pricing ker- 
nels is underway. More details will be given elsewhere ^. 



8 




0.98 I ' ' ' ' ' 1 

0.0 100.0 200.0 300.0 

days 

Figure 1. An example of thermalized path obtained from the simulation of the path 
integral (Black-Scholes) with r=0.05 and cr = 0.12 
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xO = log(SO) 



Figure 3. Call option price for strike price 3 versus the logarithm of the initial value of 
the stock XQ = log(5o). the parameters are fixed as in figs 1. Shown is the analytical 
result vs the monte carlo result, with a low resolution of 10,000 configurations 
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Figure 4. Plot of the option price versus the stock price obtained by a Langevin simu- 
lation of the path integral, with a potential step 



